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The coordinate Bethe ansatz solutions of the XXZ model for a one-dimensional spin-1/2 chain are analyzed 
with focus on the statistical properties of the constituent quasiparticles. Emphasis is given to the special cases 
known as XX, XXX , and Ising models, where considerable simplifications occur. The XXZ spectrum can 
be generated from separate pseudovacua as configurations of sets of quasiparticles with different exclusion 
statistics. These sets are complementary in the sense that the pseudovacuum of one set contains the maximum 
number of particles from the other set. The Bethe ansatz string solutions of the XXX model evolve differently 
in the planar and axial regimes. In the Ising limit they become ferromagnetic domains with integer-valued 
exclusion statistics. In the XX limit they brake apart into hard-core bosons with (effectively) fermionic 
statistics. Two sets of quasiparticles with spin 1/2 and fractional statistics are distinguished, where one set 
(spinons) generates the XXZ spectrum from the unique, critical ground state realized in the planar regime, 
and the other set (solitons) generates the same spectrum from the twofold, antiferromagnetically ordered 
ground state realized in the axial regime. In the Ising limit, the solitons become antiferromagnetic domain 
walls. 
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1. Introduction 

Quantum spin chains are physically realized in quasi-one-dimensional magnetic compounds. 
These are crystalline materials with magnetic ions arranged in exchange-coupled chains that are 
isolated from each other by non-magnetic ions. The elementary magnetic moments are localized 
on the sites I of a one-dimensional lattice, which makes them distinguishable. The associated spin 
operators thus commute if they belong to different sites of that lattice, [S^,S^,] = iheap-yS'^ 5w . 
The Hilbert space of a spin-1/2 chain with N sites is conveniently spanned by product basis vectors 
|(Ti . . .(Tat), cr/ =T,i. 

Prominent among the many models employed in the context of quantum spin chain compounds 
is the spin-1/2 XXZ model, 

N 

Hxxz = JY. {SfSf+, + 5f + A5f 5f+i } . (1) 
(=1 

It describes a uniform nearest-neighbor exchange coupling with uniaxial anisotroy. Periodic bound- 
ary conditions are assumed. We distinguish ferromagnetic coupHng (J < 0) from antiferromagnetic 
coupling (J > 0), and the planar regime (0 < A < 1) from the axial regime (A > 1). Important 
special cases are the XX model (A = 0), the XXX model (A = 1), and the Ising model (A — > oo). 

One persistent challenge through decades of experiments on quantum spin chain compounds has 
been the interpretation of the observed intensity spectrum in terms of constituent quasiparticles 
[[U [2]. The chief motivation of the work reported here is to shed new Hght on this very issue. 
The approach taken is eclectic in nature, combining older and more recent results for the unified 
purpose of understanding the quasiparticle composition of the XXZ spectrum more thoroughly. 

The strong interest in the XXZ model is sustained not only by its direct relevance in exper- 
imental physics, but also by its amenability to exact analysis via Bethe ansatz. The solution of 



© Ping Lu, Gerhard Miiller, and Michael Karbach 



1 



Quasiparticles in the XXZ model 



the XXX model was, in fact, the very problem for which Bethe originally invented the method 
in the early days of quantum mechanics [ [3l |4] . The Bethe ansatz allows for a characterization of 
all many-body eigenstates as composed of quasiparticles that scatter off each other or form bound 
states and thus turn into different quasiparticles. Over the years the same basic idea has been 
successfully appHed to many different kinds of systems and situations including interacting boson 
and fermion gases with contact interactions [lUISllT], the Hubbard model for electrons on a lattice 
with on-site repulsion [HIISI, and for bosonic and fermionic quantum field theories including the 
quantum Sine-Gordon and Thirring models [llOj. 

The XXZ spectrum can be generated from different pseudovacua by the systematic creation of 
quasiparticles with different exclusion statistics. In most cases these pseudovacua are states of low- 
est energy (physical vacua) . An overview of the different kinds of quasiparticles that have emerged 
from analytic work on the XXZ model is shown in Fig. [TJ The boxes in the top row represent 
quasiparticles whose names are derived from jargon used to describe coordinate Bethe ansatz so- 
lutions. These quasiparticles are characterized by strings of complex momenta with common real 
part and different imaginary parts [[7]. One-strings are unbound magnons, two-strings are bound 
magnon pairs etc. Small (large) imaginary parts indicate loose (tight) binding [HI]. At A = 1 the 
strings are constituent particles of (degenerate) multiplets of eigenstates with total spin St ■ These 
multiplets split up energetically at A 7^ 1. 



XX 


XXX 


Ising 


broken strings, 




strings 




stretched strings, 


fermions 






domains 



7 



spinons 




solitons 





Figure 1. Zoo of quasiparticles that play some role in the context of the XXZ model. For J > 0, 
the pseudovacua of spinons and solitons are at the bottom of the spectrum. The pseudovacuum 
of the string particles is at the top in the axial regime (A > 1) and moves downward in the 
planar regime (0 < A < 1), reaching the center in the XX limit and coinciding there with the 
pseudovacuum of Jordan- Wigner fermions. 



At A = the energy levels join up in new degenerate configurations, refiecting different sym- 
metries [ O [131 [TH [m [TBI [17] • A. more natural classification scheme for the XX spectrum is then 
based on string fragments. Hence the name broken strings. These fragments are closely related to 
the free lattice fermions that emerge from the Jordan- Wigner representation of spin-1/2 operators [ 
[T2I [T5I [m [17] . At A > 1 the imaginary parts of the string solutions grow in magnitude and diverge 
as A ^ 00. Hence the name stretched strings. The tightly bound strings are related, in the Ising 
limit, to localized domains of reversed spins on successive lattice sites. The exclusion statistics of 
these domains is similar to yet subtly different from that of the strings. 

The boxes in the bottom row of Fig. [T] represent quasiparticles that are complementary to the 
string particles. The string pseudovacuum contains the maximum number of spinons or solitons. 
The spinon and soliton vacua contain strings at maximum capacity. Whereas the string parti- 
cles have integer-valued exclusion statistics, the spinons and solitons are realizations of fractional 
statistics [[18]. They are both semions but with different pseudovacua. The exclusion principle for 
semions is, roughly speaking, halfway between those appHcable for fermions and bosons. If it takes 
1/ g particles to lower the number of orbitals in a band available for occupancy by one, then g — 1 
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describes fermions, 5 = 1/2 semions, and 5 = bosons (as a limit). 

The string particles and the semionic particles are natural building blocks for the systematic 
construction of a complete XXZ eigenbasis from different pseudovacua. The configurations of 
string particles and semionic particles are both unique and preserved in every XXZ eigenstate. 
The ground state (physical vacuum) of the XXZ antiferromagnet (J > 0) coincides with the 
pseudovacua of the semions. The spinon vacuum is unique and coincides with the non-degenerate 
ground state of the XXZ model at < A < 1 for even N. The soliton vacuum is twofold and 
coincides with the ground state of the XXZ model at A > 1 for N oo. 

Let us briefiy illustrate the relation between the particles from the top and bottom rows in Fig.[T] 
with two simple scenarios, one in configuration space and the other in momentum space. In the 
Ising limit, Hxxz has simple product eigenstates. The ferromagnetic state | ItT ' ■ ' ) is the unique 
pseudovacuum for domains of consecutive fiipped spins, H ■ ■ ■ |, as shown in top part of Fig. [21 The 
twofold Neel state, | tit ' ' ' )) I iti ' ' ' )) by contrast, is the pseudovacuum for antiferromagnetic 
domain walls of the kind IT; ii as shown in the bottom part of Fig.O These domain walls (named 
solitons) have effective spin ±1/2. The integer- valued exclusion statistics of domains is associated 
with the fact that neighboring particles must be separated by any positive integer number of lattice 
sites. The fractional exclusion statistics of domain walls, on the other hand, is associated with the 
fact that one of two lattice sites may be shared by neighboring particles. 
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Figure 2. Constituent quasiparticles of Hxxz at A — * cxa in real space (lattice of A'^ sites). They 
are either domains of flipped spins embedded in the unique vacuum with all spins up (top) or 
domain walls embedded in the twofold vacuum with all spins alternatingly up/down (bottom). 



In the XX limit, Tixxz is equivalent to a system of Np free lattice fermions in a band [ 
[Tot [201 EB [22] • The ground state corresponds to the intact Fermi sea at \k\ < kp as shown in 
Fig.[3l The entire spectrum can be generated systematically via particle excitations {ANp = +1), 
hole excitations {ANp — —1), and particle-hole excitations {ANp — 0). The Fermi-sea ground state 
can be viewed as the pseudovacuum for semionic spinous. We introduce a threshold momentum 
kc that varies with Np. When a fermion is removed from the intact sea, kc slightly increases to 
generate not just one but two vacancies in the interval < kc of the band. These two holes 
are then identified with a pair of spin-up spinous. When a fermion is added outside the intact 
Fermi sea, kc slightly decreases to produce two particles in the region > kc. They are identified 
with a pair of spin-down spinous. A fermionic particle-hole excitation leaves kc unchanged and is 
interpreted as a pair of spinous with opposite spin orientation. The entire spectrum can thus be 
described in the form of spinon configurations. 

The main goal of this work is to identify the relationship between the string particles and the 
semionic particles for Hxxz at A = 0, 1, oo in particular and to illuminate how this relationship 
can be tracked between the three points in parameter space. In Sec. [21 we set the stage for this 
investigation by reviewing the relevant Bethe-ansatz representations that best serve our purpose. 
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Figure 3. Constituent quasiparticles oiTLxxz at A = in reciprocal space (fermionic band with 
Fermi momentum kp)- They are either spinless Jordan-Wigner fermion or semionic spinons with 
spin ±1/2. 



We then discuss the quasiparticle composition of the XXZ spectrum in the XXX limit (Sec. l3.ip . 
the Ising Hmit (Sec. [121), and the XX Umit (Sec. [331) ■ 

2. Bethe ansatz 

The coordinate Bethe ansatz for Tixxz uses the magnon vacuum \F) = | ft • • • t) with energy 
Ep = JNA/4: as its reference state. The Bethe form [S| of any eigenvector in the invariant subspace 
with magnetization Mz = N/2 — r, 

\i,)= a(ni,...,n,)S- (2) 

l<ni<...<nr<W 

has coefficients of the form 

determined by r magnon momenta ki and one phase angle, 9ij — ~Oji, for each magnon pair. The 
sum V G Sr is over the permutations of the indices {1, 2, . . . , r}. The ki and % satisfy the Bethe 
ansatz equations (BAE), 



11 ' g*(fci+fcj) _|_ 1 _ 2Ae*'=3 



The energy and the wave number of an eigenstate with magnon momenta {fci, . . . , kr} are 

E — Ep 



J 



= (cos h-A), /c = ^ fc, niod(27r). (5) 



Reasons of practicality dictate the use of different transformations of the BAE in the regimes 
of axial anisotropy (A > 1), planar anisotropy (A < 1), and isotropic exchange (A — 1). It is 
convenient to introduce rapidities Zi, 

r ^ ^ Zi 

cot — tanh -— : < A = cos 7 < 1 
^, I 2 2 

cot ^ = <^ : A = 1 , (6) 



2 



coth — tan -— : A = cosh 77 > 1 
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which renders the Hmit A ^ 1 smooth from both sides. The BAE ([4]) thus transform into 

r 

N^l{z^) = 2^Ih+^^2{z,-ZJ), i=l,...,r, (7) 

where 

I arctan cot — tann — 
V 2 2 



(\)y{z) = 2 arctan — 



2 arctan coth — tan — 



< A < 1 
A = 1 
A > 1 



(8) 



The Bethe quantum numbers (BQN) li of integer or half- integer value reflect the multivaluedness 
of the logarithm used in the transformation. They are subject to restrictions that will be discussed 
case by case. The energy expression in ^ becomes 



E-Ef 
J 



^e(zi), 



e(z) 



sin^ 7 



cosh(7z) — cos 7 
2 



1 + Z2 



sinh r\ 



cos(?7z) — cosh 77 



< A < 1 
A = 1 
A > 1 



(9) 



Note that we are, effectively, deahng with a single anisotropy parameter, 77 = 77, that is real 
in one regime, imaginary in the other, and zero at the isotropy point. This parametrization is 
particularly useful for tracking the spectrum between the axial and planar regimes across the 
point of higher rotational symmetry. Slightly different parametrizations are more adequate for the 
exploration of the limits A ^ cx) in the axial regime and A ^ in the planar regime. 

2.1. Axial regime 

At A > 1 we use the transformation [[23]. 



Z'i , 1 ^ 

tan — — tanh — cot — , 
2 2 2' 



77 = arcosh A, 
instead of ([61 , to bring the BAE ([H into the form 



-TT < < TT, 



N 



coth(77/2)tan(z,,/2) - i 
coth(77/2) tan(zj/2) + 1 



The associated trigonometric BAE are 



n 



coth(?7) tan [(z.; — Zj)/2] — i 
coth(?7) tan [(z^ — Zj)/2] + z' 



7V0i(zi) = 2-Kli + ^ (?!)2(zi - Zj), 



1 r 

-L, . . . , / , 



where 



, . ^ . „ ( tan(z/2) \ 

0v\z) — 2 arctan — -— 

' Vtanh(?7;//2) / 



27r 



3?z 
27r 



(10) 



(11) 



(12) 



(13) 



The second term in (fT3|) ensures that the set remains the same as an eigenstate is tracked 
across the axial regime [[3. Here \x\ is the floor function (integer part of x). 

For the analysis of solutions that include complex magnon momenta we invoke the string 
hypothesis for the rapidities [l^: 



z^-' =z^ + 777(771+ 1-20, ^ = 1, 



771 = 1,2, 



(14) 
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The index I distinguishes rapidities belonging to the same string (of size to). The index « = !,..., rim 
distinguishes different strings of the same size. The string ansatz lfT4|) produces only asymptotic 
solutions of the BAE (fTTI) for N oo. The finite-iV corrections are, in general, exponentially 
small, and not all finite- solutions fit the string template lfT4|) [[?]. However, neither corrections 
nor exceptions affect macroscopic systems. In the Ising Hmit, where the spread of imaginary parts 
in (dH) diverges, all corrections and exceptions disappear even for finite N. 

A given string solution of ifTTj) with magnetization Mz = N/2 — r is described by r rapidities 
that breaks down into configurations of strings such that the constraint. 



^ TO rim = r, 



(15) 



is satisfied, where the set C identifies those sizes of strings that occur in a given eigenstate. With 
the functions 

^ coth(?7t//2)tan(z/2)-z _ sin ((z - ir]:^)/2^ 



co\h{^v|2) tan(z/2) + i (^^^ ^ ^^^^ ^3) 



we rewrite Eqs. (fTTj) in the form 



), « = 1, 



When we substitute the string ansatz lfT4|) into Eqs. (fTTl we obtain 



11 VlK^cc )11^2(^;a' - Z^' ), 

k^l 



(m',/3)^(m,a) k—1 

for ? = 1 , . . . , TO and a = 1 , . . . , , where 

sin (^[z™ + I'qim + l-v- 2l)]/2^ 



sin + ^??("^ + l + v- 20]/2^ 



To determine the real parts, z™, we form the product of all Eqs. ifTSl) for fixed to, a: 



L;=i 



n 

{m' ,a.' )^{rn.a.) 



/ m .A;\ 



1=1 k=il 



.,1 



(16) 



(17) 



(18) 



(19) 



(20) 



Each expression in square brackets can be simplified massively, producing the BAE for the z™ 



N 



95m'-m(2aa' )'Pm'+m(^aa' ) H [V'r.i' +m-2l(Zaa ) 



(21) 



with z:^^ 



z™ . The associated trigonometric BAE, 

7V</>,4z™)=27rC+ ^ $™n'(C"'), 

(m',Q')7^^(7n,a) 



(22) 



^ ^j-^^ ^ |0|m-m'|(^:) + 2(/!)|m_„'|+2(z) + ••• + 2(/)„i+m'-2(z) + (?!)m+„i'(z) : m' m ^^^^ 
1 2^2(2;) + 2^4(z) + ... + 2^2m-2(z) + '^2m(z) : m' = TTl ' 
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depend on a set {/™} of BQN that reflects the speciflc string combination of any given eigenstate. 
The energy and wave number of that state are 



E — E]^ 



E 

(m.a) 



sinh?7sinh(rym) 
cosh(77TO) — cos z™ 



(m.a) 



2tt 



mod(27r). 



(24) 



The range of the /™ will be discussed flrst for the case A = 1 in Sec. 12.21 and then for the axial 
regime including the Ising limit in Sec. [ 



2.2. Isotropic exchange 

At A = 1 we retain the rapidities from ([6]), the trigonometric BAE in the form and the 
energy expression ([9]). The string hypothesis now reads 



m,l rn 



z™ + i(m + l-20, 



1 = 1,. 



TO = 1,2, . 



(25) 



(26) 



The BAE for the real parts z™ are Eqs. ([ST]) with |(T9]) replaced by 

^ i ^ z^^ + i{m+l- 1^ - 21) 
'^'^^^"''^ z™ + i(TO+l + i/-20' 

The associated trigonometric BAE for the then take on the form (|22l) with ^mm' (z) from l(23|) 
and (puiz) = 2arctan(z/i/) from ([8]). The energy and wave number of an eigenstate specifled by 
the set {/™} are 



E — Ep 



E 

(7n,a) 



2to 



J m? + (z™)2^ 



k = 



(7n,a) 



mod(27r). 



(27) 



The string hypothesis sets the range of the /™ (with r^+i > ^™ implied) as follows [[7]: 
1 



E' 

ra'^C 



"mm' ' ^m' 



2min(m, m') — (5^, 



(28) 



where Um is the number of TO-strings (distinguished by running index a) in the eigenstate. Note 
that the range of allowed values becomes narrower for all sizes if a string of any size is added. The 
I™ for a given combination {rim} and a given value of to are either all integers or all half- integers 
such that the border values of the range l(28|) are realized. 

The Bethe state characterized by a set {/™} is the highest- weight component (i.e. the state 
with Mz = St) of an S'T-multiplet with total spin 



mriri 



(29) 



The other components of any given S'T-multiplet are generated by the addition of magnons with 
zero momentum. These have no effect on the energy or the wave number. Consider the case of a 
highest-weight state with total spin l(29|) that only contains 1-strings. Suppose this state is specified 
by the following set of BQN subject to the constraint (|28| : 



-l{N-r-l)<ll<---<I^<^{N 



(30) 



The rapidities z^^ derived from the BAE |[2l 



(31) 



7 



Quasiparticles in the XXZ model 



are real and, with rare exceptions, finite. The member state with = St — 1 of the same multiplet 
has one extra rapidity, z^+i — oo, representing the additional magnon with K+i — 0. The BAE 
l(22|) for this state are then satisfied with the same z^, a — 1, . . . , r and with z^+i = ±oo if we set 
the BQN as follows: 

/i=/i±i, a = l,...,r; /i+i = ±^(iV - r). (32) 

In the more general case, where the highest-weight state under consideration contains strings with 
m > 1, the shifts in the already existing BQN, /™, and the value of the new BQN, I^+i, will be 
different. In Sec. 13.11 we will treat the strings as interacting particles and examine their exclusion 
statistics. 

2.3. Planar regime 

The string ansatz at A < 1 [[7] will not be used here. Simplifications and residual complications 
that occur in the limit A ^ can be seen in the raw form Q of the BAE. There are two categories 
of solutions, both of which are ubiquitous. Regular and singular solutions are distinguished by the 
absence or presence of pairs of critical magnon momenta with ki + kj ~ tt, which make both the 
numerator and the denominator in (0]) vanish as A — > 0. 

All regular solutions produce real magnon momenta from e'^'^' — (—1)''^^, whereas singular 
solutions include critical magnon pairs that are either both real or form a complex-conjugate pair 
[ [121 [HI [TBI [E] • The latter can be interpreted as fragments of strings that exist throughout the 
planar regime. Since all critical pairs are associated with a twofold degeneracy of eigenstates with 
equal wave number, the singular features as imposed by the BAE in the limit A ^ can be 
removed by unitary transformations. This erases, at A = 0, all traces of the string nature, at 
A > 0, of complex solutions. The magnon momenta thus regularized are 

k^=7r^^ll (33) 

The associated BQN are integers for odd r and half- integers for even r with range 

\li-rr\<liN~l), r. = i[l- (-!)'■]. (34) 

Their relation to the BQN l(28|) and (|32l) of the ^T-multiplet states will be discussed in Sec. 13.31 

The regularized BAE solutions (|33l) describe hard-core bosons. The phase shift is % = tt for 
all two-particle interactions. These hard-core bosons are equivalent to the Jordan- Wigner fermions 
[ [HI [iQl [m [22] that have been instrumental in most studies of the XX model. 

3. Quasiparticle composition of XXZ spectrum 

We are now ready to explore the relationship between the complementary quasiparticle com- 
positions of the XXZ spectrum. We begin in Sec. 13. II with string particles and the complementary 
spinon particles for the XXX case. In Sec. 13. 21 we then discuss the effects of axial anisotropy on the 
strings and their relationship to ferromagnetic domains in the Ising limit. Complementary to these 
domains are the soliton particles in the shape of antiferromagnetic domain walls. In the planar 
regime the strings evolve differently. What remains of them in the XX limit are fragments that act 
like hard-core bosons or, equivalently, free Jordan- Wigner fermions. Complementary to the latter 
are again the spinous as will be discussed in Sec. 13.31 

3.1. XXX limit: strings and spinons 

The very structure of the coordinate Bethe ansatz suggests that the strings l(25|) can be inter- 
preted as quasiparticles. There exists a universal energy- momentum relation as implied by l(27l) . 
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The particle interaction is encoded in the set of momenta (or rapidities) dictated by the BAE and 
in the phase shifts associated with elastic two-particle collisions. 

The exclusion statistics of strings is determined by the rule l(28|) governing the range of BQN 
and by the relation ((29| governing the capacity for strings in a highest- weight state of given St- The 
total number of S'T-multiplets with string content {ni,n2, ■ ■ ■} becomes the solution of a standard 
combinatorial problem: 

W{{n„i})^ TT ( V dm^Am^y^ gmm'(nm' - 5mm'), (35) 

mec ^ ' m'ec 

where 

Am ^ N + 1 - 2m, gmm' — 2 min(TO, to') (36) 

are statistical capacity constants and statistical interaction coefficients, respectively, that are spe- 
cific to the string particles [ [131 [HI HH US [26] . Taking into account the {2St + l)-fold degeneracy 
of each multiplet, this classification accounts for the complete spectrum, 

^ W^({n„})(25T + l) = 2^, (37) 

{n„.} 

with the dependence of St on {rim} given in (|29l) . 

To illustrate the string composition of the XXX spectrum and to explain its relationship to 
the complementary spinon composition we consider a chain of iV = 6 sites. In Table [T] we list all 
combinations {um} permitted by (|29l) . Also listed is the range of all I™ present in each combination 
as inferred from l(28|) . The distinct configurations of /™ thus allowed produce the highest-weight 
states of all S'T-multiplets as shown in Fig. [4| (left). The quantum number St depends on the 
combinations {n^} via (|29|) and the wave number k depends on the configurations {/™} via 
l(27j) . Each row of /™ in Fig. [4] (left) represents the distinct string motif of an ST-multiplet. The 
solution of the BAE l(22|) thus specified is for the highest- weight component of the multiplet. The 
non-highest-weight components are characterized by additional BQN as described at the end of 
Sec. [221 The role of these additional BQN will be discussed in Sec. [321 for A > 1 (see Fig. [5]) and 
in Sec. [331 for A < 1 (see Fig.^. 

Table 1. String composition {rim} of all S'T-multiplets for TV = 6. Each row describes a distinct 
combination for a total of 7. Each combination produces W multiplets for a total of 20. Each 
multiplet represents 2St + 1 states for a total of 64. Associated with each string in a given 
combination is a BQN 7™ of the range shown. 
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Elements of the string motif serve as the template of the spinon motif of the same S't -multiplet. 
Spinons are specified by their spin and momentum quantum numbers. The spinon interaction 
depends on the particle momenta and (in general) also on the particle spins. The unique XXX 
ground state for even TV is the spinon vacuum. The ground state for odd N is fourfold degenerate 
and contains exactly one spinon. The total number Ng of spinon is restricted to be even (odd) for 
even (odd) N and can assume the values < TV^ < TV. Counting the spin-up spinons (iV+) and 
spin-down spinons (TV-) separately produces the relations 

TV++TV_=TV,, TV+ -TV_ = 2T\^^. (38) 



9 



Quasiparticles in the XXZ model 



-3 -2 -10 1 2 3 
I I I I I I I I LJ I LJ 



L 



® 



® 



® 



® 



® 



® 


® 






® 




® 




® 






® 




® 


® 






® 




® 






® 


® 



® ® ® 
® 



r 



2 

® 
2 

3 I 



® 



1 



-3 -2-10123 
I I I I I I I I I I LJ I 



J ^ H a a a ^ 
a H a a 



o 

a o 

a a 

a a a o 

a a a a 



a a a 
o a a 
a 

o 



o 


o 


a 


a 


-1.2808 


o 


a 


o 


a 


-1.0000 


o 


a 


a 


o 


0.1180 


a 


o 


o 


a 


-2.1180 


a 


o 


a 


o 


-1.0000 


a 


a 


o 


o 


-1.2808 


a 


a 


a 


a 


0.7808 


a 


a 


a 


a 


0.5000 


a 


a 


a 


a 


0.7808 



o a a 
a o a 
a a o 

^ a a ^ 



E/J 

1.5000 

1.0000 
0.0000 
-0.5000 
0.0000 
1.0000 



-2.8028 

-0.5000 
-1.5000 
-0.5000 

0.8028 



Figure 4. Specification of all S'T-multiplets of the XXX model for A'' = 6 in the string repre- 
sentation (left) and spinon representation (right). The quantum numbers St and k = Nk/2-K 
of each multiplet are stated on the far left. The positions of the numbers m = 1, 2, 3 and the 
vertical lines (left) mark the values of the 7™ and their range. The emphasized by circles, 
are markers for the spinon configurations (right), where they are reproduced as small circles. 
Each square on the right represents a spinon. The set of squares is complementary to that of 
circles with the same range. The number mi = 1, 2, . . . , 5 inside each square marks the spinon 
orbital. The energy of each St multiplet is stated on the far right. 



It is useful to introduce spinon orbitals associated with distinct spinon momentum quantum 
numbers m^. The wave number of any XXX multiplet can be expressed in terms of the spinon 
orbital momenta as follows: 

/„ ^= AT A 

mod(27r). (39) 




The range of rm depends on both N and Ng : 

™. = i^.i| + 2,...,iV-f. (40) 

The number of available orbits to spinons, iVorb ^ [N — Ns)/2 + 1, thus decreases by one for 
every two spinons added, AA^orb = ~^Ns/2, which is one way of recognizing the semionic nature 
of spinon particles. The exact exclusion statistics of spinons is encoded in the number of ways 
spin-up spinons and N- spin-down spinons can be distributed among iVorb accessible orbitals. The 
spinon multipHcity expression [ ._L8j , 

W{N+,N_) = n f " M ' d^^A,-J2 9aa'{N.. - 5,,,), (41) 
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_ 1 



(42) 



is cast in the same general form as its string counterpart l(35|) . Summation of W{N^,N_) over all 
allowed values of N± accounts for all 2^ states. 

The complementary relationship between the string and spinon particles is reflected in their 
motifs. In Fig. |4] we show the string motif and spinon motif of all 20 iST-multiplets for TV = 6 side 
by side. The number of spinons contained in each eigenstate of a given S'T-multiplet is equal to 
the number of vacancies left by the 1-string BQN, across the range l(28|) . In Fig. |4] (right) we 
have marked the positions of the l\ by circles and the vacancies by squares. The number inside 
each square denotes the orbital rrii to which every spinon belongs. The available orbitals depend 
on Ns via ()40| . 

The rules for assigning spinons of particular momentum quantum numbers to orbitals must 
ensure (i) that the rm reproduce, via (|39l) . the wave number k already known via l(27l) and (ii) 
that the permissible spin orientations are consistent with the quantum number St- The allowed 
combinations of spinon orbitals are encoded in those S'T-multiplets that do not contain any strings 
with m > 1. Here the positions of the 1-strings, i.e. the circles in the spinon motifs, play the role 
of delimiters between successive spinon orbitals. 

A spinon orbital with U spinons has an orbital spin S'f'^ — li/2, i = l,...,iVorb- If all Ng 
spinons are in the same orbital such as in the first two motifs (numbered from top to bottom) 
this spinon configuration represents a multiplet with St = S™^ = Ns/2. Any distribution of 
spinons into multiple orbitals represents more than one S'T-multiplet. The multiplets represented 
are determined via the decomposition of tensor products of orbital spins. 

For example, if we have one spinon with nii = 2 and three spinons with rrii = 4, that decom- 
position reads | ® | = 2 © 1. The associated multiplets with St = 2 and St = 1 are found in 
motifs three and thirteen, respectively. The configuration with two spinons in each of the same two 
orbital represents three ST-multiplets according tol(8)l = 20l0O. They are found in motifs 
four, fourteen, and twenty, respectively. The systematic application of these rules assigns a unique 
spinon momentum and spin content to every S'T-multiplet. 

The energetic split of different ST-multiplets associated with the same spinon momenta is caused 
by a coupHng between orbital spins. In the Haldane-Shastry model, which has higher symmetry, 
that coupling is absent and these particular ST-multiplets remain degenerate [ [27l [23 . In the XX 
model, which has lower symmetry, even the ST-multiplets split up energetically as will be discussed 
in Sec. [3:31 [[14]. 

The translation and refiection symmetries of Hxxz dictate that every state with wave number 
k can be transformed into a state with wave number 2i: — k and the same energy, implying that 
any state with /c = 0, tt is its own image. These transformation properties are refiected in the string 
motif (/^ — > — /™) and in the spinon motif {nii —^N — rrii) as is evident in Fig. [4l 

3.2. Ising limit: stretched strings, domains, and solitons 

In the following we describe three ways of generating the spectrum of the Ising chain, each 
producing a distinct set of quasiparticles with different exclusion statistics. We begin with the 
analysis of the string solutions of the BAE from Sec. 12.11 To arrive at a non-divergent spectrum in 
the Ising limit of Hxxz we rescale the exchange coupling as follows: 



JV 




(43) 



The rescaled energy expression ^ of the Bethe ansatz in raw form thus becomes 




(44) 
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Evidently, real magnon momenta contribute only summarily to the energy, namely via the first 
term in (|44l) . Non- vanishing terms cos ki/ A can only come from complex ki with infinite imaginary 
parts. 

We first track the highest- weight states of the S'T-multiplets from the XXX limit to the Ising 
Hmit. For these states, all solutions that are real at A = 1 stay real and all solutions that are 
complex at A = 1 have imaginary parts that diverge as A ^ oo. The non-highest-weight states 
evolve far less uniformly. The additional rapidities, which are all equal to ±00 at A = 1, become, 
in general, finite and thus contribute to the energetic split of the ^r-multiplets. In some states, 
the extra Zi stay real, in other states they combine to form complex pairs. It appears that this 
transformation from real to complex rapidities takes place throughout the axial regime. Here we 
focus on the end product at A = cxo, where the ferromagnetic domains in Ising product eigenstates 
become the natural quasiparticles. 

Expanding the trigonometric BAE (|22|) of the axial regime about the Ising limit produces, in 
leading order, a set of linear BAE for the real parts of the rapidities. 




2nJ^- ^ (45) 



m'ec 



where the new set {J™} of BQN, related to the original set {/™} by a shift that depends on N 
and Um, has the same range l(28|) 



|-^n<^(^^l^Xl ) ' = 2min(TO,r7i') - Smm' , (46) 

\ m'GC / 

which guarantees the correct number of states from the highest-weight type. As the first step in 
the solution of the Hnear BAE l(45|) for these states we introduce the quantities 



m' a— 1 

Equations (|45l) summed over a can thus be brought into the form 



27r- 2^ -— imm'^m' (48) 

^™ m'ec 



with C,m guaranteed to be positive. Equation (|48l) is solved by matrix inversion. Substitution of the 
solution Zm into l(45|) yields the rapidities 

^ (49) 

m'ec 

More expHcit solutions in compact form for the situations where only 1-strings or only 2-strings 
are present read 

1 27r , 27r ^ 
Zr. = — ■Ja~TuTT Q; = l,...,r, (50) 



z 



2 



7V-r " N{N-r) 
27r ^9 67r 



, Ji T^J2, a=l,...,r/2. (51) 

iV-3r/2 " Ar(iV-3r/2) ' ' ' ' ^ ' 

If a single m-string with 1 < m < A^/2 is present we have z™ = (27r/7V) J™ with | Jf | < {N — 2m)/2. 

In Fig. [H we have reproduced the BQN of all 20 highest-weight states for TV = 6 at A = 1 from 
Fig. [4] and have added the BQN for all 22 non-highest-weight states with Mz > 0. The additional 

of non-highest- weight states, obtained for A = 1 as described at the end of Sec. 12.21 are located 
beyond the range of the for the highest-weight states (vertical lines) . 

The choice of the extra and the associated shift of the /™ already present is not unique in 
most cases. Our choice was guided by the aim to avoid multiple of the same value and to restore 
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Figure 5. Specification of all XXX eigenstates with Mz > for A*' = 6 in the string representation 
(two columns). The quantum numbers M^, St and k = Nk/2TT of each state are stated on the 
left. The vertical lines indicate the range l[28ll of the /„. The positions of the numbers m = 1, 2, 3 
between the vertical lines represent the values of the J™ for the highest-weight states (with 
= St)- The additional BQN of the non-highest-weight states (with < Mz < St) are all 
located on the outside. In the Ising limit each state evolves into a translationally invariant linear 
combination of Ising product states with domains of lengths as indicated on the far right in each 
column. 

the symmetry relations described at the end of Sec. 13.11 For states with A: = 0, tt, which are their 
own images under the symmetry transformation, neither goal was fully achievable for iV = 6. The 

configurations are asymmetric in several instances and two identical are unavoidable in one 
instance. These aesthetic fiaws pose no problem for what follows. 

To the right of each string motif in Fig. [5] we have added another motif in the shape of a ket. It 
identifies each string solution as tracked from A = 1 to A = oo with an Ising eigenstate as composed 
of Hnear combinations of ferromagnetic domains. We have already introduced these domains in the 
context of Fig. [21 Here we use a notation conducive to linking Ising product eigenstates with obvious 
domain composition to Bethe eigenstates with specific string structure in the limit A = oo. 

Let the set of A product vectors with r — N/2 — fiipped spins that are generated from one 
such state, |(Ti • ■ • CAf), via translations be represented by {|(Ti • ■ • (Tm)}\- Let 

{|1_^2_^3_^---)}^ (52) 

Ul V2 V3 

represent the set of A translationally invariant linear combinations of all product states that contain 
Vn domains of length n subject to the constraint z/i -I- 2:^2 + 3:^3 -I- • ■ ■ — r. Any set of the latter kind is 
constructed from one or several sets of the former kind with matching domain content. For = 6 
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there are seven sets of translationally invariant Ising eigenstates at > thus constructed: 

{|0)}i {| TTTTTT)}i; {|i)}6 {| TTTTTDIe; {|2)}6 {| TTTTiDle; (53a) 
{|ii)}9'^ {| TTTiTi>}6,{| TTiTTDIa; {|3)}6 {| TTTiiDle; (53b) 
{|i2)}i2 {\ TTiitDle, {| TTiTii)}6; {|iii)}2 {\ TiTin)}2. (53c) 

In Fig. [5] we have identified each of these 42 states with a particular BAE solution in the 
limit A — > oo. The identification is straightforward for the highest-weight states. Each BQN /™ 
represents exactly one domain of size = m. However, the evolution of the non-highest-weight 
states between A = 1 and A ^ oo is far less predictable. The extra BQN are all of the type 
at A = 1 and the associated rapidities are real and infinite. In some cases all extra rapidities stay 
real. The domain structure of these states is determined by the /™ just as in highest-weight states. 

However, for the majority of non-highest-weight states some of the rapidities that start out 
real at A = 1 turn into complex-conjugate pairs at A > 1 with imaginary parts that diverge as 
A ^ oo. This metamorphosis necessitates a reconfiguration of the associated BQN (not shown in 
Fig. [5]). Of the 22 non-highest- weight states for iV = 6 only six do not acquire additional complex 
rapidities between A = 1 and A oo, namely those states for which the string content encoded 
in the /™ matches the domain content encoded in the ket. 

The identification of the BAE solutions in the limit A ^ oo with Ising eigenstates of specific 
domain content may be less certain in longer chains where states with different {i^^} but equal Mz 
are degenerate. That is the case for the two sets of states {|13)}24 and {|22)}i2 in a chain with 
N = 8 sites, for example. However, it seems reasonable to assume that this degeneracy is absent 
at A < cx) in most if not all cases and thus guarantees that the BAE solutions have, in general, a 
unique domain content in the limit A ^ oo. Even though the systematics of the transformation of 
some BAE solutions between A = 1 and A — > oo appears elusive at present, the systematics of the 
endproduct, namely the domain structure of the Ising spectrum for arbitrary N, is well known. 

An Ising chain of length N can accommodate domains with = 1, . . . , N — 1 consecutive fiipped 
spins. Domains of size fj, are treated as distinct species of independent particles. The capacity for 
domains is restricted by the relation 

N-l 

^ (/i + 1)1.^ < N. (54) 

The total number of Ising eigenstates with domain content {i^i, i'2, . . .} is governed by a multiplicity 
expression somewhat similar to l(35|) yet different [[26]: 

= n " ~ ' d,^A,-J2 5mm'(^m' - S,,'), (55) 

Ai=i ^ A" / 

A, = N-,, 9,,^[ l;^^^ , r^j:^... (56) 

The capacity of the system for domains is controlled by and the statistical interaction between 
domains by g^^'. All domains have the same energy, J. The statistical mechanics of domains for 
the Ising chain has been carried out exactly [|2Q, reproducing familiar results. 

In the context of Fig. [2] we have already qualitatively described the antiferromagnetic do- 
main walls (solitons) that are complementary to the ferromagnetic domains. Among the four dis- 
tinct bonds in the general product state \(Tia2 ■ ■ ■ <y tq) , the bonds IT, ii represent solitons with 
spin -1-1/2, —1/2, respectively, and ti, it are vacuum bonds. Close-packed solitons with like 
spin orientation reside on successive bonds (e.g. TTT); whereas close-packed solitons with oppo- 
site spin orientation are separated by one vacuum bond (e.g. ttii)- Each of the seven sets of 
Ising eigenstates l(53|) contains a specific number of spin-up and spin-down solitons, (A^+,iV_) = 
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(6, 0), (4, 0), (3, 1), (2, 0), (2, 2), (1, 1), (0, 0). All solitons have the same energy, J/2. The multiplicity 
expression for solitons [I26j. 

^(^+'^-) = ^^r^ n { '^"^^^'''^ V d, = A,- Y.g„.,{N,.-5..,), (57) 

A„^]^(N-l), w = Na = N++N^, (58) 

is similar to expression (|4T1) for spinons. Solitons and spinons are both semions but have different 
pseudovacua. Whereas the spinon vacuum was found to be unique, the soliton vacuum is twofold, 
consisting of the two product Neel states, | Tit ' ' ' i) and | iti ' ' ' T) or linear combinations 
thereof. The soliton vacuum, like the spinon vacuum, is realized only in chains with even N . 

3.3. XX limit: broken strings, fermions, and spinons 

Here we investigate what happens to the highest-weight and non-highest-weight states of the 
^T-multiplets in the presence of planar exchange anisotropy, particularly in the Hmit A = 0. 
We have seen that in the axial regime there exists a tendency for real BAE solutions to become 
complex. As A ^ cx) all imaginary parts diverge, binding the magnons tightly into domains. In the 
planar regime, there exists a trend in the opposite direction. At A = all BAE solutions can be 
regularized as stated in Sec. 12.31 making all magnon momenta real. All strings with m > 1 break 
up into 1-strings. The configurations allowed by the set l(34|) of BQN, which describe regularized 
solutions, produce the complete spectrum for arbitrary TV. 

In Fig. [6] we show the configurations of of all 42 eigenstates with > at A = for 

= 6 in the same sequence as the corresponding S'T-multiplet states at A = 1 have been listed 
in Fig. [5l All multiplet states at A = 1 that contain only 1-strings keep the same configuration 
at A = 0. In the other states, all stay in the same position except the pair of that have 
identical values in Fig. [H They are replaced in Fig. [6] by a pair of distinct Furthermore, each 

in Fig. [5] is replaced by two l\ in Fig. [6] and the one occurrence of a is replaced by three 

Whereas the general rules for these substitutions still elude us, the constraints imposed by 
symmetry and conservation laws eliminate any ambiguity for the case iV = 6 shown here. The 
resulting configuration of produces, via (|33l) . the exact configurations of Jordan- Wigner fermion 
momenta in a periodic chain [IHlEl]. The statistical mechanics of the regularized 1-strings is that 
of a system of free Jordan- Wigner fermions [ [J^ . 

Each XX eigenstate thus identified by its composition of 1-strings has a unique composition 
of spinons. The rules for inferring the spinon motif from the string motif are consistent with the 
rules described in Sec. l2.2l for the case A = 1 but there are some noteworthy differences. The lower 
rotational symmetry at A = 0, which splits up the S'T-multiplet degeneracy, makes it possible to 
assign to each XX state not only a spinon configuration with unique momentum quantum numbers 
but to the spinons in each orbital also a unique spin orientation. 

The spinon motif is encoded in the string motif of each state shown in Fig. [6] as described in 
the following, (i) Consider the vertical lines dividing the space of the into two domains the 
inside and the outside, the latter wrapping around at ±N/2. (ii) Every /^-vacancy (small circle) 
inside represents a spin-up spinon (marked by a square) and every ('1') outside represents a 
spin-down spinon (marked by a diamond), (iii) Any number of adjacent spinons in the motif are in 
the same orbital, i.e. have the same momentum quantum number rrii from the set l(40|) . Two squares 
or diamonds that are separated by i consecutive 'I's have spinon quantum numbers separated by 
2£. (iv) The spinon momentum quantum numbers are sorted in increasing order from the line on 
the left toward the right through the inside domain {mf ) and toward the left with wrap-around 
through the outside domain (m~). 

The spinon motif of Fig. [6] establishes the much needed link between spinon motif of the 
fermion representation at A = introduced in Refs. [[131 E] and the spinon motif of the string 
representation at A = 1 introduced in Fig. [5l Unlike the solitons at A — > cx), the spinons at A = 
are not free. Nevertheless, an exact statistical mechanical analysis of spinons at A = is possible 
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Figure 6. Specification of all XX eigenstates with Mz > for A*' = 6 in the string representation 
(two columns). The quantum numbers Mz, and k = Nk /2n of each state are stated on the left. 
Also shown on the left is the quantum number St of the multiplet at A = 1 from which the 
state evolved. The states are listed in the same sequence as those in Fig. O The positions of 
the 'r represent the values of the 1^ from l|34p and constitute the motifs of the (broken) string 
configuration of the XX eigenstates. The motif of the complementary spinon configuration is 
marked by the squares (spin-up spinous) and diamonds (spin-down spinous) as explained in the 
text. The momentum quantum number rrii and the spin orientation o"; of each spinon present in 
that state are given to the right of the dual string/spinon motif. 



[ fT4]. Away from the Ising limit, the soUtons become interacting particles as well. They can be 
tracked all the way from A^cx)toA = O.A comparison between the spinon composition and 
soliton composition of the XX eigenstates can be found in Ref. [ [26] . 

4. Conclusion 

The long-established integrability of Hxxz imposes stringent constraints on the nature of 
the interaction between quasiparticles of any kind, reducing it, effectively, to two-body scattering 
events. As a consequence, these quasiparticles have infinite lifetimes and can thus be regarded as 
structural elements of the XXZ eigenstates. This makes it possible to systematically generate the 
complete spectrum from the pseudovacuum of one or the other set of quasiparticles. 

In this work we have discussed two complementary sets of quasiparticles, one set being the basic 
elements of string solutions of the coordinate Bethe ansatz and the other the semionic spinous or 
solitons. Our focus has been on developing interlocking motifs for each set of quasiparticles for 
the purpose of tracking them across the axial and planar regimes of Hxxz- Starting from the 
symmetry point A = 1 in parameter space we have identified opposite trends, in the two regimes, 
in the evolution of string particles, which can be thought of as bound clusters of magnons. In 
the axial regime, the strings show a tendency of increasing tightness in the binding within each 
cluster (rapidities with growing imaginary parts) and a tendency of merging clusters (real rapidities 
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becoming complex). The endproduct at A = oo are ferromagnetic domains of flipped spins. In the 
planar regime, the strings appear to loosen up and break apart. At A = it is possible to transform 
away all complex rapidities and with them any trace of magnon clustering. The free magnons behave 
like hard-core bosons or free fermions. 

For J > the pseudovacuum of string particles is located at the top of the spectrum throughout 
the axial regime and then migrates down toward the center of the spectrum in the planar regime. 
The ground state (physical vacuum) of Tixxz is unique in the planar regime and twofold in the 
axial regime. It can be identified as the pseudovacuum of spinons or solitons, respectively. We have 
given detailed descriptions of spinons at A = and A = 1 and of solitons at A = oo, including 
their relationship to the string particles. What remains to be investigated is the metamorphosis 
of the string particles between the special points A = 0,1,cxd in parameter space. This is work 
in progress. Of particular interest is an analysis of the transformations for which we have given 
here only qualitative descriptions and the impact of these transformations on the motifs such as 
displayed in Figs. |4l[6] including the effects on the complementary spinon and soliton particles. 
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